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By using the partial transpose and realignment method,we study the time evolution of the bound 
entanglement under the bilinear-biquadratic Hamiltonian. For the initial Horodecki's bound en- 
tangled state, it keeps bound entangled for some time, while for the initial bound entangled states 
constructed from the unextendable product basis, they become free once the time evolution begins. 
The time evolution provides a new way to construct bound entangled states, and also gives a method 
to free bound entanglement. 
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Quantum information theory has drawn much atten- 
tion in the last decade, and in the quantum information 
theory, quantum entanglement plays an important role 
as it can be utilized to realize some quantum processes 
In realistic world, due to the interaction between 
the system and the environment, most states in nature 
are mixed and the set of mixed states is dense in Hilbert 
space. The mixed entangled states can be divided into 
two classes 0], one is free, which means that the state 
can be distilled, the other is bound, which means that the 
state cannot be distilled. The bound entangled state is 
thought to be useless in quantum communication. How- 
ever, the bound entangled state (BES) can produce a 
nonclassical effect, which is called activation of bound 
entanglement. The underlying concept originates from 
a formal entanglement-energy analogy 0, [j, 0] • It im- 
plies that the bound entanglement is like the energy of a 
system confined in a shallow potential well. If we add a 
small amount of extra energy to the system, its energy 
can be deliberated. One of the main consequence of the 
existence of BESs is that it reveals a transparent form 
of irreversibility in entanglement processing Q . The ir- 
reversibility can be viewed as an analog to irreversible 
thermodynamics processes [7]. 

The bound entanglement can possibly get free if we add 
some energy to it. In this paper, we see what happens 
if we let the BES undergoes some quantum dynamical 
evolution under physical Hamiltonians. We will see that 
BESs can be free via time evolution. The time-evolved 
BES can be also considered as the original BES with 
an extra time parameter t, and this provides a way to 
construct more general BESs. 

There are some measures to detect entangled state 
in high-dimensional Hilbert space. For example, Peres- 
Horodecki criterion based on partial transpose (PT) 
[1, [13] and the realignment criterion [ll|, [l^] ■ Consider 
a density matrix p in a n<S)n system , the PT with respect 
to the second system, and the realignment of the density 
matrix is given by [p^'^)rj,ki = P%i,k] [p^)i].ki = Pik,ji , 
respectively. Two quantities are defined as 
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where the first is the negativity [10]. The trace norm 
II A \\ is given by ]] A ]]=trVAAt. Either TVi > or 
N2 > indicates that the state is entangled, A^i = and 
> indicates that the state is bound entaiigled, and 
A^i > means that the state is free entangled flS] . 

To study time evolution, we should have a suitable 
Hamiltonian. We will concentrate on 3 ® 3 system in 
this study, and without loss of generality, regard this 
system as a two spin-one system. The most representa- 
tive interacting model in spin-one system is the bilinear- 
biquadratic model [l4| . The corresponding Hamiltonian 
is given by . 



= Si • S2 - /3(si • S2)" 



(2) 



Here, denotes the spin-1 operator. As we know, for 
spin-1 system, the swapping operator S12 can be written 

as fii 



S12 = (Si • 82)^ + Si • S2 - 1, 



(3) 



which is invariant under the SU{3) unitary transforma- 
tion. The singlet projection operator is written as 



Pl2 = ^[(si-S2f -1] 



(4) 



Then, the Hamiltonian can be rewritten in terms of swap- 
ping operator and projective operator. 



i/ = Si2-3(l-l-/3)Pi2-/3. 



(5) 



It is obvious that swapping operator and projective 
operator are commutative, so that the unitary matrix 
U (t) can be obtained as 



U{t) = e"*-^* g-itSl2gi3(l + /3)tPl2g»i/3^ 

We have 8^2 = 1 and PJ2 = P12, then 

g-itsi2 ^ cosi - isiniSi2, 
gj3(l-(-/3)tPi2 ^ 1 + (e*3(i+'3)* - l)Pi2- 



(6) 



(7) 



(1) 



In this way, we can get the evolution operator exactly. 
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For clarity, we consider the case of /? 
tonian can be reduced to the swap 



H = S = S 
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-1. TheHamil- 



(8) 



up to a additive constant. In the following, we will study 
the time-evolution of bound entangled state govenered 
by the Hamiltonian. The initial state will evolve under 
the unitary density matrix 



U{t) 



-itH 



cost — isiniS. 



(9) 



At time t, the density matrix p(t) 



-itH 



p(0)e 



itH 



The two-spin BES we are considering, described by 
Horodecki in Ref. 0], is given by 



Pa(0) 

where 
P 



-£»-, 2 < a < 5, 



(10) 



Q+ 
Q- 



= |*+)(*+|,|M/+) = — (|00) + |11) 
= i(|01)(01| + |12)(12|-f |20)(20|), 



|22)), 



-(|10)(10| + |21>(21| 



|02)(02|) = Sg+S. 



In Ref. Horodecki have demonstrated that 

separable for 2 < a < 3 
Pa is \ bound entangled for 3 < a < 4, 
free entangled for a < 4 < 5. 

An interesting feature of this state is that it is invariant 
under the joint transformation of exchange a <-*■ 5 — a 
and swap of two particles. It is clear that the density 
matrix pa has four distinct eigenvalues, which are 2/7, 
0, (5 — a)/21, q;/21. They are non-degenerate, 2-fold 
degenerate, 3-fold degenerate, and 3-fold degenerate, re- 
spectively. 

From the expression of the state, after PT, the trace 
norm of the partially transposed state is given by 



1 



2+ -(5+ ^41 - 4a(5 - a) 



5- ^/4^^^AcJ^^l) ) 
1 for 2 < a < 4, 



i 2+ ^41-4a(5 



for a < 4 < 5. 



Note that the realignment in the realignment criteria 
is equivalent to a one-side swap followed by a PT p^ . 



namely, 



Spa = yP+ + ^Sg+ 
Then, after PT, we have 



It is direct to show that 
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g+S. 



^ (l9 + 2v/l9 - 3a(5 - a)) . 



From the expression of the two trace norms, an observa- 
tion is that the exchange a < — > (5 — a) does not change 
the values of the two trace norms. 

Now we consider the time evolution of state pa ■ From 
the evolution operator ^ and the initial state pU|) . the 
state at time t is obtained as 



Pa{t) = i{2P+ + a(<)^7+ + 6(^)^?_+^c(^)[S,^)_]X,ll) 

a(t) — a cos^ t + {5 ~ a) sin^ t 
b(t) — a sin^ t + {5 — a) cos^ t 

c(t) = sinicosi(2a-5) = l/2sin(2t)(2Q-5) (12) 

In the derivation of the above equation, we have used the 
identity [S, p_] = [g+, S]. It is obvious that 



, + 6 = 5, a - 6 = (2q; - 5) cos(2i). 



(13) 



which will be used later. 

Let's study the case in which t = 7r/2. At this moment. 



U{t = n/2) = 
So the density matrix becomes 



(14) 



PaiTr/2) = U{Tr/2)paU^{Tr/2) 
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a 
7 



g~ 



7 



-g+- 



(15) 



The above state can be obtained from Pct(O) by just ex- 
change a < — > (5 — a). But the exchange does not affect 
the values of trace norms. Thus, the trace norms at time 
t = tt/2 are the same as those at time t = 0. 

The negativity can be calculated by first making the 
PT, and diagonalizing the partially transposed density 
matrix. Making use of the following properties. 



iP-Sf^ 
we have 



3' 
1 

3 
1 



g± = g±, 



(|00>(11| + |11)(22| + |22)(00|) 
(|11)(00| + |22)(11| + |00)(22|), (16) 



Pi%t) 



1 

21 



{2S + 3a{t)g^ 



-zc(t)(|00)(ll| 
-zc(t)(|ll)(00| 



+ 36(t)£»- 

|11)(22|- 
|22)(11|- 



|22)(00|) (17) 
|00)(22|)}.(18) 



In the basis spanned by 

{|00), 111), |22), |01), |10), |12), |21), |20), |02>}, matrix 
Pa^ i^) '^^^ be written in a block-diagonal form 
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By. 



2 ic(t) -ic{t) 
= I -ic{t) 2 ic{t) 
ic{t) -ic{t) 2 



,B2 



a{t) 2 
2 6(t) 



where we have used the fact that in this basis, the swap 
can be written as 

S = /3x3 ®(rx®(^x® err- 
and /g^g is a 3x3 identity matrix, and cix is the x com- 
ponent of the Pauh matrix vector. 

Having the block-diagonal form of p-^^^^^^gne can ob- 
tain the eigenvalues of the partial transposed density ma- 
trix, and only the two following eigenvalues are possibly 
negative 



Ai 



= f 2 - %/3c 



21 



^(4-V3| sin(2t)(2a-5)| 



A, = 



-(a(t) + h{t) - 
1 



— (^5 - v/16 + (2a - 5)2 cos2 2tj , (19) 



So, the negativity of Pa{t) is given by 

Ni = max(0, -Ai) -f 3 max(0, -A2). 



(20) 



Having known analytical expression of A'^i, we next 
consider the quantity A^2- First, we need to obtain matrix 
[Spait))^''. From the definition of the EES, we have 

Spa(0 = ]^{2V+ + a{t)Sg+ + b{t)Sg-+ic{t) {q_ - q+)}. 

After PT, we obtain 
1 



{Spo.{t)) 



T2 



2^{2S + ^3c(^)(^7_-^?+) 

+a{t) (|11)(00| + |22)(11| + |00)(22|) 
+6(t)(|00)(llK|ll)(22| + |22)(00|)}, 



where Eq. is used. In the basis given above, the 

matrix {S pa{t))'^^ can be written as 

{Spa,it)f^ - ^ (^3,3 e i?2x2 © ^^2X2 © ^^2x2) , 



2 b{t) a{t) , ^ _■ (,\ 



Again, it is of block-diagonal form. The square roots of 
eigenvalues of CC^ are given by 




A/l + 3(2a-5)2cos2(2i), 



(21) 
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FIG. 1: Time evolution of the negativity A^i and the quantity 
versus t with initial state being Horodecki's state. 



and the square roots of eigenvalues of DD'^ are 



Ci = |2 + c| = -|4 + sin(2f)(2a-5)|. 



C2 = |2 



-|4-sin(2i)(2a-5)|. 



(22) 



Thus, one has the trace norm 



J + ^Vl + 3(2a-5)2cos2(2t) 



1 
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|4-|-sin(2i)(2a - 5)| 



|4-sin(2<)(2a-5)|). 



(23) 



From the expression of the trace norm, when 2 < a < 4.5, 
the above expression can be simplified as 



\P 



^ (19 + v/l + 3(2a-5)2cos2(2t)) 



(24) 



for any time. Quantity N2 is obtained by substituting 
Eq. dMl) to (P). 

Now, we analyze the analytical results of negativity 
iVi and quantity 7V2. One can see that when t — 7r/2, Ni 
and N2 are both back to the value when i = 0, because 
the swap operator has no effect on the entanglement of a 
bipartite system. When t = tt, the state has undergone 
a cycle. One can find when a — 5/2, all the Ai keep 
constant with the time. This is due to the fact that the 
BES with a = 5/2 is invariant under swapping operation. 

From Eqs. and (gD]), it is not difficult to find that, 
for arbitrary time, A^i = when 2 < a < | + --2=« 3.655. 

Then, we know that in this range if N2 > 0, the state is 
bound entangled. In Fig. 1, we plot the negativity A^i 
and the quantity N2 versus t for different a. When the 
initial state is a separable state {a = 2.8), one cannot 
generate entanglement. For a = 3.5, the initial state is 
bound entangled, and as time increases, the strength of 
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FIG. 2: Time evolution of the negativity A''i and the quantity FIG. 3: Time evolution of the negativity Ni and the quantity 
N2 versus t with initial state being first BES constructed from N2 versus t with initial state being second BES constructed 
UPB. from UPB. 



bound entanglement decreases until it vanishes. After 
some time, the state becomes bound entangled again. 
Here, there is no free entanglement. 

When we increase a to a — 3.9 > 3.655, the state 
can be made free entangled over a certain range of time 
within one period. After time evolution begins. The en- 
tanglement keep bounded until it becomes free. From 
the figure, we also see that when the initial state is free 
entangled (a = 4.5), the entanglement keeps free all the 
time. We see that the bound entanglement can be free 
by the time evolution. 

Now, let us see the evolution of concurrence. We can- 
not get concurrence, but can know a lower bound. In 
Ref. [131 , Chen et al. derived a theorem that for any 
mixed quantum state p, the concurrence C(p) satisfies 



C(p)> 



n{n — 1) 



[max (II 



Rip) 



(25) 



The theorem gives the lower bound of concurrence of 
an arbitrary mixed state. In the case of n = 3. From 
Fig. 1, we can read the time behaviors of the lower bound. 
For a = 3.9,4.5, the bound is always larger than zero, 
and it displays some singularities due to the competition 
between the two trace norms. 

Next, we consider BES constructed from the unextend- 
able product basis (UPB). The first BES from UPB is 
given by [3, 



1 

1 

1 



|o)(|o)-|i)), 

|2)(|1)-|2)), 103) 



1 
1 



(|0)-|1))|2), 



(|1»-|2))|0), 



(|0) + |1> + |2))(|0) + |1) + |2)), 



(26) 



1 4 

-(/9x9-El'^^)('^J-|)' 



(27) 



here /gxg is the 9x9 identity matrix. The second BES 
from UPB is given by Eq. ^ with ^ 



■j) = If^j) «) |w2jmod5), j = 0, ...4, 

^ ^ (27rj75),sin(2^j75),\/l + V5/2 



V 



V5 + \/5 



cos 



from which the density matrix could be expressed as 



The negativity of both these two states are zero, but 
the trace norm is given by 1.087 and 1.098, respectively, 
indicating that these two states are bound entangled. 

The time evolution of the two BESs are plotted in 
Fig. 2 and 3, respectively. There is no qualitative dif- 
ference between the two plots. The bound entanglement 
becomes free once triggering the time evolution, namely, 
the bound entanglement is very fragile comparing with 
that in the Horodecki's state. Also, we see that the lower 
bound of concurrence is determined by quantity N2 for 
all time. We can free bound entanglement by the time 
evolution. 

In conclusion, we have studied the time-evolution of 
BESs by two different operational approaches, namely, 
the partial transpose and realignment approach. The en- 
tanglement properties of time-evolution of different kind 
of initial BESs under the bilinear-biquadratic Hamilto- 
nian. The analytical results of trace norms are obtained 
for the initial state being Horodecki's BES at any time. 

For the initial Horodecki's BES, the state keeps bound 
entangled in the beginning of time evolution, while for 
the initial BESs constructed from UPB, the BES becomes 
free once the time evolution begins. The time behaviors 
are very different for different 'class' of BESs. Some BESs 
are relatively stable against time evolution, and some are 
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fragile. Even a small perturbation can free the bound 
entanglement from UPB. The time evolution provide a 
new way to construct BES, and also gives a method to 
free bound entanglement. 
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